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Abstract 

Current publicly available computer programs calculate the spectrum and couplings of the minimal supersymmetric 
standard model under the assumption of TJ-parity conservation. Here, we describe an extension to the SOFTSUSY pro- 
gram which includes /^-parity violating effects. The user provides a theoretical boundary condition upon the high- 
scale supersymmetry breaking 7? -parity violating couplings. Successful radiative electroweak symmetry breaking, 
electroweak and CKM matrix data are used as weak-scale boundary conditions. The renormalisation group equations 
are solved numerically between the weak scale and a high energy scale using a nested iterative algorithm. This paper 
serves as a manual to the /^-parity violating mode of the program, detailing the approximations and conventions used. 
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1. Program Summary 

Program title: SOFTSUSY 

Program obtainable from: http : //proj ects .hepf orge . org/ sof tsusy/| 

Distribution format: tar.gz 
Programming language: C++, f ortran 
Computer: Personal computer 
Operating system: Tested on Linux 4.x 
Word size: 32 bits 
External routines: None 

Typical running time: a few seconds per parameter point. 

Nature of problem: Calculating supersymmetric particle spectrum and mixing parameters in the 7?-parity violating 
minimal supersymmetric standard model. The solution to the renormalisation group equations must be consistent with 
a high-scale boundary condition on supersymmetry breaking parameters and Ijip parameters, as well as a weak-scale 
boundary condition on gauge couplings, Yukawa couplings and the Higgs potential parameters. 
Solution method: Nested iterative algorithm. 

Restrictions: SOFTSUSY will provide a solution only in the perturbative regime and it assumes that all couplings of 
the MSSM are real (i.e. CP-conserving). The iterative SOFTSUSY algorithm will not converge if parameters are too 
close to a boundary of successful electroweak symmetry breaking, but a warning flag will alert the user to this fact. 
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2. Introduction 



Spectrum generators are a widely used tool in particle physics beyond the Standard Model (SM), especially in 
the case of Supersymmetric (SUSY) models. Spectrum generators can be used in theoretical studies of a SUSY 
breaking scheme, for example in studies of fine-tuning. Phenomenological investigations of new patterns of SUSY 
breaking require a calculation of the spectrum. Often, the resulting SUSY spectrum is used to calculate the prospects 
of high energy experiments such as the Large Hadron Collider (LHC) discovering and measuring SUSY particles, 
assuming some SUSY breaking scheme fT] . In order to run a realistic collider and detector simulation of a new physics 
signal, a consistent model is needed as input. Such simulations are required in order to set search and measurement 
strategies as well as to estimate SUSY backgrounds to some measurement. In the event of discovery of some 
SUSY signals in LHC data, attention will turn to the question of which patterns of SUSY breaking are consistent with 
data. In such tests, SUSY spectrum generation would be an essential step. SUSY studies often perform parameter 
scans, resulting in a large number of generated spectra. There is thus a need for accurate and quick computation of 
the supersymmetric spectrum as a first step in such studies. There exist several publicly available spectrum generators 
for the /^-parity (Rp) conserving minimal supersymmetric extension of the Standard Model (MSSM): ISASUGRA 13], 
SOFTSUSY |5], SUSPECT 1^ and SPHEND |7]. Spectrum information is typically transfered to decay packages and 
event generators via a file in the SUSY Les Houches Accord format |8, 9]. 

The most general renormalisable superpotential of the MSSM contains /^-Parity violating (ftp) couplings, violat- 
ing baryon and lepton number BlOtl . A symmetry can be imposed upon the model in order to maintain stability of the 
proton, for example baryon triality fl)^ or proton hexality 1 12]. It has been shown that^p models may have interesting 



features, such as the generation of non-zero neutrino masses without the addition of right-handed neutrino fields ||13|] . 
and the gravitino as a viable dark matter candidate il4ll . The violation of baryon or lepton number implied by Ijtp leads 
to additional possibilities for SUSY detection, since such quantum numbers are conserved in the perturbative SM. 
There are important implications for direct collider searches, since one can lose the classic large missing transverse 
energy "smoking-gun" signature of SUSY? All of these features make I^p worthy of study. There is thus a strong 
motivation to extend the /?p-conserving spectrum generating public computer codes to include l^p effects. Here, we 
describe such an extension which has been applied to SOFTSUSY. The latest version of SOFTSUSY including effects 
can be downloaded from address 



http : //pro j ects . hepf orge . org/ sof tsusy/| 



Installation instructions and more detailed technical documentation of the code may also be found there. 

The Rp conserving aspects of SOFTSUSY are already explained in detail in ref. |5], and so they shall not be 
repeated here or throughout this manual, which will concentrate solely on the I^p aspects of the calculation. Adding 
ftp couplings roughly doubles the (already large) number of couplings of the MSSM. The calculation is thus more 
complicated and so it takes considerably longer than the Rp case (roughly a factor of three for identical precision). 
However, this still means that a single point in parameter space can be calculated in a couple of seconds on a modern 
personal computer The added complication of means that some features of the ftp version of SOFTSUSY are less 
accurate the the /J^-conserving case, using only one-loop RGBs to evolve the couplings and masses of MSSM fields, 
as opposed to two-loop RGBs in the Rp case. Therefore, taking the^^ calculation in the limit of small couplings, the 
numerical values of SOFTSUSY outputs will not exactly agree with the 7?p-conserving version of SOFTSUSY. We stress 
though, that if the user does not desire to include l^p couplings, the code automatically uses the Rp calculation with 
the associated speed and accuracy. Where the accuracy of SOFTSUSY in the ^p-mode differs from the /?p-conserving 
mode calculation, we shall make a note. 

We proceed with a definition of the SOFTSUSY convention for the I^p parameters and mixings in section[3] Next, 
in section m we discuss the calculation, making a note of parts which differ in accuracy to the SOFTSUSY Rp cal- 
culation. Installation instructions can be found on the SOFTSUSY web-site, but instructions to run the program can 



be found in Appendix Appendix A The output from a SOFTSUSY sample run is displayed and discussed in Ap- 



pendix Appendix B whereas a sample main program is shown and explained in Appendix Appendix C Some more 



technical information on the structure of the program can be found in Appendix [Appendix D It is expected that the 



information in Appendix Appendix D will only be of use to users who wish to 'hack' SOFTSUSY in some fashion. 
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3. MSSM Parameters 



In this section, we introduce the Ijip MSSM parameters in the SOFTSUSY conventions. The translations to the 



actual variable names that are bein g us ed in the program code are shown explicitly in appendix Appendix D The Ijip 
SOFTSUSY calculation follows ref. lll5ll and so the notation and conventions are similar. 



3.1. Supersymmetric parameters 

The chiral superfield particle content of the MSSM has the following S U(3)c'xS U(2)iXU{1)y quantum numbers: 



L: (1,2,-^), E: (1,1,1), e:(3,2,^), f7:(3,l,-^), 
D: (3,1, i), Hr. (1-2,-i), //2: (1,2,1). 



(1) 



L, Q, Hi, and H2 are the left handed doublet lepton and quark superfields and the two Higgs doublets. E, U, and 
D are the lepton, up-type quark and down-type quark right-handed superfield singlets, respectively. Note that the 
lepton doublet superfields and the Higgs doublet superfield coupling to the down-type quarks. Hi , have the same 
SM gauge quantum numbers. The l^p part of the renormalisable MSSM superpotential is written, in the interaction 
eigenbasis. 



Wi 



(2) 



Here, we denote an SU{3) colour index of the fundamental representation by e {1,2,3). The SU(2)l fun- 

damental representation indices are denoted by {a,b,c] e {1,2) and the generation indices by {i,j,k] e {1,2,3). 
e"^'^ and eab - e"'' are totally antisymmetric tensors, with £123 = 1 and ei2 = 1, respectively. Currently, only real 



^xyz 



couplings in the superpotential and Lagrangian are included. 



3.2. ftp SUSY breaking parameters 

We now detail the notation of the soft.^i'p SUSY breaking parameters. The trihnear.^Jp scalar interaction potential 



IS 



y. 



- £ah 



1 



-hijtLlL''jh + h'.j,L^Qfdk, + H.c. 



^^xyzh^Xd^dl + H.c. 



(3) 



where fields with a tilde are the scalar components of the superfield with the identical capital letter The electric 
charges of u, d, and e are -|, |, and 1, respectively. "H.c." denotes the Hermitian conjugate of the preceding terms. 
The bilinear scalar interaction potential is given by 



m? „ V: m + H.c. 



(4) 



5.5. Tree-level masses 

The mixing of MSSM particles can change in the case that lepton number is violated by the l^p interactions. Two 
cases of lepton number violating mixings are implemented in SOFTSUSY: neutrino-neutralino mixing and chargino- 
lepton mixing. We neglect sneutrino-anti-sneutrino mixing, because this has been shown to have negligible phe- 
nomenological consequences once experimental bounds have been applied 1^6*1 . 

In the presence of lepton number violating l^p interactions, the neutrinos mix with the neutralinos. At tree level, 
this results in one massive neutrino, two massless neutrinos and four massive neutralinos. ^^-loop corrections to the 
neutral fermion mass matrix (currently neglected by SOFTSUSY) can result in all neutrinos acquiring masses and the 
emergence of a PMNS mixing matrix in lepton charged current interactions. The (7 x 7) neutrino-neutraUno mass 
matrix for the three generations of neutrinos is given in OlSll and reads 



I? 
"2 



(5) 



where 
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(6) 



where /<•,■ are the biUnear mixing parameters in Eq. (|2]i, v, are the sneutrino vacuum expectation values (VEVs) and 
the Ml , M2 are the gaugino masses of hypercharge and weak isospin, respectively. The matrix ^ has five non- 
zero eigenvalues, i.e. four neutralinos and one neutrino. We denote the mass eigenstates which are obtained upon 



diagonalisation of M (in ascending order of mass): v,=i 2,3,jt'i 234^ with masses along the diagonal of the matrix 



(7) 



where (9 is a member of (9(7). A simple multiplication of rows of O by factors of / can absorb any minus signs in 

In addition, charged leptons mix with the charginos. The Lagrangian contains the (5 x 5) lepton-chargino mass 
matrix 



£^ -(-i^-,hl,e, ,,)Mc 



+ H.C. 



(8) 



The mass eigenstates £ - (e,/i, t),Xi2 ^re given upon the diagonalisation of the matrix 

0; 



Mc 



( M, ^v„ 

13.,,. 



(9) 



Here, Ye is the lepton Yukawa matrix from the Rp superpotential in ref . isll ■ We define the diagonalised mass matrix 



(10) 



U and V being orthogonal 5 by 5 matrices. 



4. Calculation Algorithm 

In broad terms, the algorithm for the calculation of the Ijip MSSM spectrum follows that of the 7?^ -conserving 
case, although some of the individual steps differ in the Ijtp case. It is performed via the iterative algorithm depicted 
in Fig. [1] An initial estimate of gauge couplings and up quark masses at Mz are obtained as in the Rp version of 
SOFTSUSY in ref. |5]. However, the charged lepton and down-quark Yukawa couplings receive corrections, and are 
detailed in section|4T| The MSSM parameters are then run to the scale 

Ms = ^|mr,{Ms)n%(Ms\ (11) 

where the scale dependence of the electroweak breaking conditions is small |[l3l. Electroweak symmetry breaking 
(EWSB) conditions are then imposed, as described in section 14.31 taking into account sneutrino VEVs and other 
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1. SUSY radiative corrections to gi(Mz) 



2. Ye, Yd include v^^ contributions. Iterative solution of Ye- 



convergence 



3. Run to Ms. 



r 



4. EWSB, iterative solution of fi and sneutrino VEVs. 



5. Run to Mx- Apply soft breaking and^,, SUSY boundary conditions. 



6. Run to Mz- 



L 



7. Calculate Higgs and sparticle pole masses. Run to Mz- 



Figure 1: Iterative algorithm used to calculate the Ijip MSSM spectrum. The initial step is the uppermost one. Ms is the scale at which the EWSB 
conditions are imposed, as discussed in the text. Mx is the scale at which the high energy SUSY breaking boundary conditions are imposed. 

lepton-number violating effects. The MSSM parameters are then run up to some high energy scale Mx, where the soft 
SUSY breaking terms are fixed by a user-specified boundary condition. Mx may be the electroweak gauge unification 
scale, some scale pre-specified by the user or indeed Mj. The running of the MSSM couplings is described in 
section l4!2l The SUSY Ijtp couplings Ajjk, /l^^, /tj^^, are then fixed at this scale Mx- The model is then run down 
to Mz, when the iteration is performed again by returning to step 1 in Fig. [T] Iteration proceeds until, at step 3, all 
parameters evaluated at Ms are identical to within a fractional accuracy of TOLERANCE to the previous iteration's (at 
step 3). TOLERANCE< 1 is a numerical parameter set by the user, with default value 10"^. Once this has been achieved, 
the algorithm proceeds to step 7, where the pole masses of sparticles are calculated as in section l4~4l 

4.1. Gauge and yukawa couplings 

In order to calculate the Yukawa couplings of the down quarks, contributions to the mass matrix from sneutrino 
VEVs are taken into account: 



where all parameters are evaluated at Mz and are in the DR scheme in the MSSM. The down-quark mass matrix in 
the weak eigenbasis, {mD)ij, is obtained as in the Rp SOFTSUSY version. Vd is the VEV of H\, as obtained below. 

The chargino-lepton mixing in Eq. (|9]l, complicates the matching of {Ye),] to the charged lepton masses. We 
employ an iterative procedure in order to calculate which {Ye))] predict the empirical input values of charged lepton 
masses. 

1. Initially, we set (if),; as in the -conserving limit, ignoring any charged lepton-chargino mixing or sneutrino 
VEVs. Thus, the empirical MSSM DR values of charged lepton masses evaluated at Mz are written in the 
matrix (m''^'')ij as diagonal values. Then, (Ye)ii - {m'^'')u '^/vd (no sum on 1). 

2. The resulting matrix (Fe),-; is then to substituted into Eq. (|9]l to form the charged lepton-chargino mass matrix, 
obtaining the 5 x5 U and V transformation matrices that diagonalise it via Eq. (fTOl i. 




(12) 
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We denote the 3 by 3 lower right-hand side blocks of U and V hy U and V, respectively. The Yukawa matrix 
Ye is then set to be 

(13) 



At present, physical lepton mixing is not implemented in the public version of SOFTSUSY. 

This result for Ye is then inserted back into step 2, leading to a better approximation of U and V. Steps 2 to 

4 are iterated until successive iterations predict identical diagonal entries of Ye within a fractional accuracy of 

lO^'^XTDLERANCE. 



4.2. Running ofMSSM couplings 

For the TJ^-conserving parameters, the renormalisation group evolution (RGB) employs two-loop MSSM p func- 
tions for the supersymmetric parameters (l^, including tan/3 and the Higgs VEV parameter v. Gaugino masses and 
/Jp-conserving SUSY breaking Higgs parameters are also run to two-loop order in the i?p-conserving parameters. 
The other Tip- conserving SUSY breaking parameters (sfermion mass matrices and some tri-linear couplings) may 
be set to two-loop or one-loop order by the boolean parameter INCLUDE_2_L00P_SCALAR_CDRRECTI0NS in the main 
program. The RGB includes full family dependence and the complete set of 1-loop MSSM I^p /3 functions in both 
SUSY-preserving and SUSY-breaking l^p parameters 01511 . The increased number of ftp couplings and yS- functions, 
as well as other complications, mean that the I^p mode runs more slowly than the 7?^ -conserving mode. In the case 
of running SOFTSUSY in the f(p mode the accuracy does not match the extremely high one of the Rp version, in order 
to keep the running time down. In either the -conserving or the Ijtp mode, the program can be made to run faster 
by switching off the two-loop renormalisation of the scalar masses and tri-linear scalar couplings. All /3 functions are 
real and include inter-generational quark mixing effects. 



4.3. Electroweak symmetry breaking 

We now discuss the minimisation the potential of the neutral scalar fields, {/I'i, Zip vi, V2, V3) at the renormalisation 



scale Ms- Following the calculation in Ref. UlSll . this system of equations is solved using the following definitions 

Vu 



tan 



n ^" 2 2 2 2 



The VBVs Vd and v„ can be written 



(14) 



= cos2/^(v^-|]vf), v^ = sin^/3(v^-gv?). 



(15) 



We see from Bq. ( fTsT i that the presence of sneutrino VBVs does not change the numerical value of tan^S. This 
convenient formulation was first developed in Ref. [|19il . The BWSB condition for the Higgs superpotential mass 
term yu can be written 1 15] as 



1 



tan2/3 - 1 



—2 / 2 X ^1' . '^i 



1 



2n.,2 n 



<. + kr-T(r + ^2K-A- 



tan^/J-iM^j, (16) 



The soft SUSY breaking mass squared term m^ is expressed in terms of /i, v„, Vd, and v,-: 



sin2yS 



+ m]j^ + 2\fi\^ + 



{m~ ) + k'h 



Vj ~ Vi 

— -Di- 

Vd Vu 



(17) 



where we employ the simplifying notation 



-2 _ 2 ^^d\Y 

ma = mij H — , 

Vu dvu 



—2 



Vd dvd ' 



(18) 
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The tadpoles |^ currently only contain the 7?^ -conserving contributions. Ijip contributions to them are currently 

neglected, for these are of order so they are small for small couplings Aijk,A'.jf, as is impUed for most 

{/, j, k} by the experimental bounds lEol] . In the 7?^ -conserving limit /c,, v,, D,, ^ ) ^ 0, Eqs. ( fTSI ) and ( fTTI i tend 
to the usual /J^-conserving MSSM Higgs potential minimisation conditions. 

The potential minimisation conditions for the sneutrino VEVs may be written as 



i 

where 



5AV] 

Vd + AV„ --g^}^ (19) 



(Ml)ij = (mi^)ji + + ^-Ml cos 2/? 6ij + ^^^^ sin^ (v' - - vj) 6ij. (20) 

We now detail the iterative procedure by which SOFTSUSY obtains parameters describing a minimum of the po- 
tential with the correct properties, i.e. satisfying Eqs. (fT6b.(fT7li and ( fT9] l. All of the running parameters discussed are 
evaluated at a renormalisation scale Q - Ms . 

1. For given value of tan/J, Eq. (fT4l i provides an initial estimate for v„ and in the 7?^ -conserving limit of v, = 0. 
Eqs. ( fTSI l and (fTTb are also first solved in the R,, limit, i.e. v,- - 0, /c,- = Z), = (m^ - ) = 0. This provides initial 
values for /i and niy 

2. The sneutrino VEVs v, are now obtained from the left hand side of Eq. (fT9] l by using v„, v^,// and as 
previously derived in the iterative procedure in the right hand side of the equation. 

3. The corrected values of v„, are then computed including the non-zero sneutrino VEVs v, via Eq. ( fTSl l. 

4. fi and are then obtained from the left hand sides of Eqs. (fT6b and (fTTb . The code returns to step 2 and steps 
2-4 are iterated until {v,, fi, m^) all change by less than a fractional accuracy of TOLERANCEx 10"** on successive 
iterations. 



4.4. MSSM spectrum 

Neutralino masses are calculated at tree-level as in Eq. (|6]l. All TJ^-conserving one-loop threshold corrections 
are then added. The neutrino masses are calculated by diagonalising this mass matrix, and taking the lightest three 
eigenvalues, whereas neutralino masses are defined to be the largest four eigenvalues. When the chargino masses 
are calculated by the iterative procedure in section WA\ one-loop /?p-conserving corrections to chargino masses are 
added to the two by two top left-hand corner of Eq. (|9]). All DR quantities in the mass matrices are taken at the renor- 
malisation scale Ms in the tree-level mass matrices. All other masses are calculated according to the /?p-conserving 
SOFTSUSY calculation, i.e. including the one-loop 7?p-conserving threshold contributions. 



4.5. Physics applications 

The Ijtp aspects of prototype versions of SOFTSUSYhave akeady proved useful for various studies, for example: in 
determining the different possibilities for the lightest supersymmetric particle in the mSUGRA framework 12 1[ 1221] . in 
defining benchmark points for future comparative collider studies Ii23i] . for studying neutrino mass textures ll24ll and 
for investigating Ijip mechanisms of neutrinoless double beta decay ll25i 12611 . 
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Appendix A. Running SOFTSUSY 



SOFTSUSY produces an executable called sof tpoint . x. For the calculation of the spectrum of single points in 
parameter space, we recommend the SUSY Les Houches Accord 2 (SLHA2) |9] input/output option. The user must 
provide a file (e.g. the example file included in the SOFTSUSY distribution rpvHouchesInput), that specifies the 
model dependent input parameters. The code may then be run with 

. /sof tpoint .X leshouches < rpvHouchesInput 

For the SLHA2 input option, the output will also be given in SLHA2 format. The example file provided calculates the 
same point as the mSUGRA point we give as an example below. Such output can be used for input into other programs 
which subscribe to the accord, such as PYTHIA |27] (for simulating sparticle production and decays at colliders), for 
example. For further details on the necessary format of the input file, see ref. [9J. Note, that SOFTSUSY does not yet 
support the (optional) setting of the bilinear sneutrino VEVs, these are instead fixed by Eq.[T9l It supports the setting 
of all other SLHA2 input blocks associated with non-complex 

For a quick examination of a single point in mSUGRA parameter space, the command 

. /sof tpoint . X sugra <mO> <ml2> <aO> <tanb> <mx> <sgnMu> RPVcoupling <i> <j> <k> <value> 

can be utilised. Bracketed entries should be replaced by their numerical values, where all massive parameters (otq, 
Ml/2, ^0, Mcut) should be quoted in GeV. RPVcoupling € {lambda, lambdaP, or lambdaPP) = {A, A', A") and 
{i, j ,k)G {1,2,3) specify a single non-zero GUT-scale coupling, mx denotes the scale at which the high-energy 
boundary condition is to be applied. If this is specified as unified, as in the Rp version, the electroweak gauge 
unification scale Mcut is used (defined to be the DR scale Q at which gi(Q) = g2(Q))- 



Appendix B. Sample Output 

For the recommended SLHA2 option, the conventions for the output are fully explained in Ref. f^. We therefore 
present the non-SLHA2 comphant SOFTSUSY output for a single example calculation. The following is an example 
for one non-zero coupling A'^22(^gut) = 0.1 at Mcut in an^^ mSUGRA SPSla-Uke 11281] scenario: 

./sof tpoint. X sugra 100 250 -100 10 unified 1 lambdaP 3 2 2 0.1 

The output obtained from this command is: 

SOFTSUSY SUGRA calculation 
mgut = 2.0970e+16 



neutrino masses : (1 ,3) : 

4.2082e-63 -1.3249e-51 -4.5924e-03 
sneutrino Vevs : (1,3) : 
-1.5184e-22 -1.0550e-16 -2 . 2648e+00 
full 7 by 7 neutral fermion mixing(7,7) : 
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OOOOe+00 
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1114e-06 
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6933e 


-23 
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0557e- 


-25 


-2 


4594e 


-25 


2 


0278e 


-26 


-4 


8200e- 


-26 


-1 


1114e-06 


1 


OOOOe+00 


4 


6474e 


-17 


2 


1172e- 


-19 


-1 


6997e 


-19 


1 


2082e 


-20 


-3 


5320e- 


-20 


-1 


5284e-23 


-4 


6475e-17 


9 


9998e 


-01 


4 


5877e- 


-03 


-3 


7151e 


-03 


4 


0463e 


-04 


-6 


2576e- 


-04 


-7 


3953e-42 


-5 


1024e-35 


-4 


1501e 


-03 


9 


8441e- 


-01 


1 


1217e 


-01 


-6 


1485e 


-02 


-1 


2067e- 


-01 


-1 


0388e-40 


-1 


3487e-35 


3 


9421e 


-03 


-6 


1991e- 


-02 


9 


4040e 


-01 


9 


1902e 


-02 


3 


2147e- 


-01 


-9 


5829e-28 


9 


6737e-21 


-1 


6184e 


-03 


1 


5415e- 


-01 


-2 


7820e 


-01 


6 


9496e 


-01 


6 


4488e- 


-01 


-8 


3721e-42 


1 


3205e-35 


1 


4448e 


-04 


-5 


7576e- 


-02 


1 


6019e 


-01 


7 


1049e 


-01 


-6 


8280e- 


-01 



RPV soft breaking parameters at Q=9.1188e+01 
Dr(l,3) : 

-3.0513e-21 5.1915e-20 -8.5363e+01 
mHllsq(l,3) : 

1.3521e-20 -2.3005e-19 3.3979e+02 
hur : 

Matrix 1(3,3) : 
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Matrix 2(3,3) 










OOOOe+00 

V/ • WWW ' \J\J 
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OOOOe+00 

, V,/ WV/W ' vv/ 
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OOOOe+00 


OOOOe+00 
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0. 


OOOOg+00 
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0. 


OOOOe+00 


. OOOOe+00 


0. 


, OOOOe+00 
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. OOOOe+00 


Matrix 3(3,3) 










. OOOOe+00 


, 


, OOOOe+00 


, 


. OOOOe+00 


. OOOOe+00 


0. 


.OOOOe+00 


0. 


, OOOOe+00 


. OOOOe+00 


0, 


. OOOOe+00 


0. 


. OOOOe+00 


hdr : 










Matrix 1(3,3) 










2.2127e-27 


3, 


, 1687e-65 


-2 , 


, 9449e-54 


-3.7647e-26 - 


■5. 


, 3912e-64 


5, 


. 0105e-53 


7.6908e-05 


7, 


,6273e-43 


-1, 


.0143e-31 


Matrix 2(3,3) 










1.0704e-63 - 


■ 1 . 


,4253e-20 


3, 


, 2565e-41 


-1 .8212e-62 


2, 


.4251e-19 


-5, 


. 5406e-40 


3 3460e-41 - 


■2 . 


8535e+02 


1 , 


0060e-18 


Matrix 3(3,3) 










-9.0631e-52 


9, 


. 5898e-40 


2, 


. 0221e-24 


1.5420e-50 - 


■ 1 . 


6316e-38 


-3. 


. 4404e-23 


-3.1480e-29 


2, 


.2915e-17 


6, 


.9662e-02 


her : 










Matrix 1(3,3) 










O.OOOOe+00 


5, 


, 1077e-27 


-9, 


,2676e-06 


-5.1077e-27 


0, 


.OOOOe+00 


-3, 


.9853e-19 


9 . 2676e-06 


3. 


.9853e-19 


0. 


.OOOOe+00 


Matrix 2(3,3) 










. OOOOe+00 


6, 


.2073e-26 


-2, 


.7304e-21 


-6.2073e-26 


0, 


.OOOOe+00 


-1, 


.9162e-03 


2.7304e-21 


1. 


.9162e-03 


0. 


.OOOOe+00 


Matrix 3(3,3) 










0. OOOOe+00 - 


■3. 


, 1239e-43 


5, 


,6720e-22 


3.1239e-43 


0, 


.OOOOe+00 


-6, 


.7578e-24 


-5.6720e-22 


6, 


.7578e-24 


0, 


.OOOOe+00 



RPV supersymmetric parameters at Q=9.1188e+01 
kappa(l,3) : 

1.7868e-24 -3.0401e-23 
Lambda_U : 
Matrix 1(3,3) : 

0. OOOOe+00 0. OOOOe+00 

O.OOOOe+00 O.OOOOe+00 

0. OOOOe+00 O.OOOOe+00 
Matrix 2(3,3) : 

O.OOOOe+00 O.OOOOe+00 

O.OOOOe+00 O.OOOOe+00 

O.OOOOe+00 0. OOOOe+00 
Matrix 3(3,3) : 

O.OOOOe+00 O.OOOOe+00 

O.OOOOe+00 O.OOOOe+00 

O.OOOOe+00 O.OOOOe+00 
Lambda_D : 
Matrix 1(3,3) : 



8.2516e-02 



O.OOOOe+00 
O.OOOOe+00 
O.OOOOe+00 

O.OOOOe+00 
O.OOOOe+00 
O.OOOOe+00 

O.OOOOe+00 
O.OOOOe+00 
O.OOOOe+00 



-8 1351e-31 - 


■ 1. 


,2976e- 


-68 


1, 


. 1127e- 


-57 


1 3841e-29 


2. 


. 2077e- 


-67 


-1 , 


. 8931e- 


-56 


-3 7569e-08 - 


3. 


, 7243e- 


-46 


5 , 


. 1385e- 


-35 


Matrix 2(3,3) 














-4 1117e-67 




. 0943e- 


-23 


-1 , 


. 2987e- 


-44 




1 , 


,8618e- 


-22 


2. 


.2097e- 


-43 


-1 .6498e-44 


3 , 


. 1408e- 


-01 


-5 , 


. 1807e- 


-22 


Matrix 3 ('3 3") 














3 3387p-B'i - 


■4. 


,0763e- 


-43 


-7 , 


.6352e- 


-28 


-5 6805e-54 


6 , 


. 9354e- 


-42 


1 , 


.2991e- 


-26 


1.5419e-32 - 


1 . 


, 1700e- 


-20 


-3. 


.5260e- 


-05 


LcinibcicL E i 














Matrix 1(3,3) 














r\ nnnriizi+nn - 
u . uuuueTVju 


z . 


,5171e- 


-30 


c 
D , 


,8320e- 


-09 


2.5171e-30 


0. 


•OOOOe+OO 


2. 


,9379e- 


-22 


-6.8320e-09 - 


■2. 


, 9379e- 


-22 


0, 


.OOOOe+00 


Matrix 2(3,3) 














O.OOOOe+00 - 


■3. 


,0589e- 


-29 


2, 


,0129e- 


-24 


3.0589e-29 


0, 


, OOOOe+00 


1, 


.4126e- 


-06 


-2.0129e-24 - 


1. 


,4126e- 


-06 


0, 


. OOOOe+00 


Matrix 3(3,3) 














O.OOOOe+00 


1, 


,5394e- 


-46 


-4, 


, 1861e- 


-25 


-1.5394e-46 


0, 


. OOOOe+00 


1, 


.3041e- 


-26 


4.1861e-25 - 


•1. 


,3041e- 


-26 


0. 


.OOOOe+00 



Gravitino mass M3/2: O.OOOOe+00 
Msusy: 4.6733e+02 MW: 8.0429e+01 
Data set : 

mU: 1.3901e-03 mC: 6.3535e-01 mt: 1.6390e+02 mfpole: 1.7140e+02 

mD: 2.7641e-03 mS: 6.0519e-02 mB: 2.8844e+00 inb(mb) : 4.2000e+00 

mE: 5.0267e-04 mM: 1.0394e-01 mT: 1.7516e+00 mb'pole: 4.9390e+00 

aE: 1.2792e+02 aS: 1 . 1760e-01 Q: 9.1188e+01 mT-pole: 1.7770e+00 



loops : 3 



thresholds : 1 



SUSY breaking MSSM parameters at Q: 9.1188e+01 



UA(3,3) : 










-5.9460e-03 


5 


. lOlOe-44 


-7, 


. 1094e-28 


1.1163e-46 


-2, 


. 6880e+00 


3, 


,9461e-19 


-6.1653e-33 


1 


.5775e-21 


-4, 


.9799e+02 


DA(3,3) : 










-1.4915e-01 


-4, 


.8416e-38 


6, 


.0862e-26 


6.6163e-45 


-3, 


.0969e+00 


1, 


,9011e-19 


1.8853e-28 


-1, 


.3852e-15 


-1. 


,2963e+02 


EA(3,3) : 










-7.3611e-03 


-2, 


. 1687e-18 


4, 


.4632e-19 


-3.1654e-16 


-1, 


.5220e+00 


6, 


.6103e-20 


4.0105e-16 


-1, 


.5223e-16 


-2. 


.3693e+01 


mQLsq(3,3) : 










3 . 5007e+05 


-1, 


. 8445e-38 


9. 


,6546e-25 


-1.8445e-38 


3 


.3981e+05 


-5, 


.3981e-16 


9 . 6546e-25 


-5, 


.3981e-16 


2. 


.9086e+05 


mURsq(3,3) : 










3.3056e+05 


7, 


.7940e-47 


-1. 


.0757e-30 


7 . 7940e-47 


3, 


.3056e+05 


2. 


.7160e-19 


-1.0757e-30 


2 


.7160e-19 


2, 


. 1378e+05 



mDRsq(3,3) : 
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3.2841e+05 -1.7248e-39 
-1.7248e-39 3.0788e+05 
8.6490e-27 -4.6844e-14 

mLLsq(3,3) : 

3.8907e+04 -2.9050e-18 
-2.9050e-18 3.8905e+04 



6.3968e-18 
mSEsq(3,3) : 
1.8572e+04 
-5.0373e-16 



-1.0883e-16 



8.6490e-27 
-4.6844e-14 
3.2398e+05 

6.3968e-18 
-1.0883e-16 
7.7821e+03 



-5.0373e-16 1.1532e-14 

1.8569e+04 -4.3773e-15 

1.1532e-14 -4.3773e-15 1.7891e+04 

m3sq: 2.4235e+04 mHlsq: 3.1426e+04 mH2sq: - 
Gaugino masses(l,3) : 

9.8542e+01 1.8879e+02 6.3121e+02 
Supersjnmnetric parameters at Q: 9.1188e+01 
Y-U(3,3) : 

1.3185e-48 -2.3776e-37 

3.6188e-03 1.1559e-23 
3.6339e-26 



1.5499e+05 



7.9176e-06 
2.8638e-51 

-3.4340e-42 
Y-D(3,3) : 
1 . 5743e-04 
2 . 0057e-49 

-2.1533e-31 
Y-E(3,3) : 
2 . 8630e-05 
1.2312e-18 

-1.7387e-18 



5.2005e-41 
3.4470e-03 
1.6089e-18 



9.0099e-01 

-6.4612e-29 
6.8201e-24 
1.4776e-01 



8.4350e-21 -1.7520e-21 
5.9198e-03 1.6806e-23 
6.5998e-19 1.0271e-01 
higgs VEV: 2.4954e+02 tan beta: l.OOOOe+01 smu: 3.4577e+02 
gl: 4.5952e-01 g2: 6.4262e-01 g3: 1.1428e+00 
thresholds: 3 #loops: 2 



Physical MSSM parameters : 

mh'O: 1.0945e+02 mA'O: 3.9171e+02 mH-Q: 

alpha: -1.1432e-01 

sneutrinos (1 ,3) : 

1.8588e+02 



3.9210e+02 mH'+-: 4.0017e+02 



1.8588e+02 
mU~(2,3) : 
5 . 6375e+02 
5.4707e+02 
mD~(2,3) : 
5.6918e+02 
5.4686e+02 
inE~(2,3) : 
2 . 0207e+02 
1.4389e+02 



5.5679e+02 
5.4707e+02 

5.6230e+02 
5.3236e+02 



1.0592e+02 

5.8667e+02 
4.0051e+02 

5.1470e+02 
5.4674e+02 



2.0208e+02 
1.4388e+02 
thetat: -5.9454e-01 thetab: 
mGluino : 6 . 0998e+02 
chargirL0s(l,2) : 

1.7934e+02 3.7564e+02 
thetaL: -4.2097e-01 thetaR: 
neutralinos (1,4): 

9.7246e+01 1.7975e+02 -3.5740e+02 
neutralino mixing matrix (4,4) : 

9.8482e-01 1.0923e-01 -6.1480e-02 
-5.8733e-02 9.3919e-01 9.1496e-02 
1.5314e-01 -2.8150e-01 6.9449e-01 



1.0975e+02 

1.5956e+02 

3.2174e-01 thetatau: 6.8797e-01 



-2.4342e-01 



3.7559e+02 

-1.2011e-01 
3.2573e-01 
6.4420e-01 
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-5.6870e-02 1.6354e-01 7.1101e-01 -6.8153e-01 
Isp is neutralino of mass 9.7246e+01 GeV 



For brevity, the numerical precision of the output has been set to 3 significant figures here, but the default option is 
for 5 significant figures. After printing the value of Mqut in GeV, the running f(p MSSM parameters are output at 
Q = Mz- The neutrino masses niy. and sneutrino VEVs v,- are output in units of GeV. Next, the 7x7 neutralino-fermion- 
mixing matrix O from Eq. (|7]l is specified. Then follow the soft breaking and l^p supersymmetric parameters from 
section [3] In order, these are: the bi-linear scalar interaction potential parameters, Z), and (m^ ), both quoted in 

units of GeV^, and the tri-linear baryon-number violating scalar interactions in terms of three anti-symmetric matrices 
(hur,), where each entry is in units of GeV. The lepton-number violating scalar interactions are also written in units 
of GeV in terms of three matrices, where 

(huri)jk = h'.'j^, (hdri)jk = h'ji^., (her^jk = hjki (B.l) 

The three /<•, parameters come next, listed in units of GeV. Finally the tri-linear superpotential interactions are listed in 
terms of the dimensionless matrices 

(Lambda_U,)jvi = A'.'j^, (LambdaJ),)^^ = .1^^,., (Lambda_E,)>i = ■Iju- (B.2) 

The rest of the information printed is standard 7?^, -conserving MSSM information from SDFTSUSY. The soft SUSY 
breaking parameters were defined in |5]. Any problems with the parameter point are printed in the final line (none in 
this case, indicating a viable point). 

Appendix C. Sample Program 

In this section we present a sample main program, that illustrates a scan over an ^p parameter This main program 
is included in the rpvmain.cpp file with the standard SDFTSUSY distribution and performs a scan in A'^^i^Mqut), 
assuming the SPSla |28] mSUGRA parameters mq = 100 GeV, M1/2 = 250 GeV, Aq = -100 GeV, tan^S = 10 and 
fi > 0. The size of the coupling varies from A'^^-^ e [0, 0.144]. The program prints out the value of the tau sneutrino 
pole mass and any problems associated with the point in question for each value of A'^^^iMcur)- 

The sample program has the following form: 

#include <rpvinain.h> 

/* global variable declaration */ 

int MIXING = -1, PRINTOUT = 0; 

bool gaugeUnif ication = true, ewsbBCscale = false; 
/// fractional accuracy required 
double TOLERANCE = l.Oe-3; 
/// decay constant of muon 
double GMU = 1.16637e-5; 

/// there are two possible conventions: if QEWSB > MZ, its value is assumed 
/// in GeV and used as a constant MSUSY. Otherwise, it MULTIPLIES the usual 
/// MSUSY value, of root(mstopl mstop2) 
double qEWSB = 1.0; 

/// Do we include 2-loop RGEs of *all* scalar masses and A-terms, or only the 
/// scalar mass Higgs parameters? (Other quantities all 2-loop anyway) : the 
/// default in SDFTSUSY 3 is to include all 2-loop terms, except for RPV, 
/// which is already slow and calculated to less accuracy than the R-parity 
/// conserving version 

bool INCLUDE_2_L00P_SCALAR_C0RRECTI0NS = false; 

/// number of loops used to calculate Higgs mass and tadpoles. They should be 
/// identical to each other for a consistent calculation 
int numHiggsMassLoops = 2, numRewsbLoops = 2; 
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/// you can define your own MZ value in GeV: this one is the default listed in 
/// the file def.h 
double MZ = MZCENT; 

/* end of global variable declaration */ 

int mainO { 

/// Sets format of output: 6 decimal places 
outputCharacteristics(6) ; 

/// Header 

cerr « "SOFTSUSY" « SOFTSUSY_ VERSION 

« " Ben Allanach, Markus Bernhardt 2009\n" ; 
cerr « "If you use SOFTSUSY, please refer to B.C. Allanach, "; 

cerr « "Comput. Phys. Commun. 143 (2002) 305, \protect\vrule widthOpt\protect\href-[http://arxiv.org/abs/hep-ph 
cerr « "For RPV aspects, B.C. Allanach and M.A. Bernhardt, "; 

cerr « "\protect\vrule widthOpt\protect\href {http : //arxiv . org/abs/0903 . 1805}{arXiv : 0903 . 1805> . \n\n" ; 

/// "try" catches errors in main program and prints them out 
try { 

/// Contains default quark and lepton masses and gauge coupling 
/// information 

QedQcd oneset; ///< See "lowe.h" for default parameter definitions 

oneset .toMzO ; ///< Rims SM fermion masses to MZ 

/// Print out the Standard Model data being used, as well as quark mixing 

/// assumption and the numerical accuracy of the solution 

cerr « "Low energy data in SOFTSUSY: MIXIWG=" « MIXING « " TOLERANCE=" 

« TOLERANCE « endl « oneset « endl; 
/// set peirameters 
double tanb = 10. ; 
int sgnMu = 1 ; 
double mgutGuess = 2.el6; 

double aO = -100.0, ml2 = 250.0, mO = 100.0; 

/// number of points for scan 
const int numPoints = 20; 

/// parameter region 

double Start =0. , End = 0.144; 

DoubleVector pars (3); 

/// set basic entries in pars 

pars(l) = mO; pars (2) = ml2; pars (3) = aO; 

cout « "1'_{331}(M_X) m_snutau # Problem flag" « endl; 

/// loop over parameter space region 
int ii; for (ii=0; ii<=numPoints ; ii++){ 

double lambda = Start + ((End - Start) / double (niunPoints) * double (ii)); 

/// define rpvSoftsusy object 
RpvSoftsusy kw; 

/// set lambda coupling at mgut 
kw.setLamPrimeO, 3, 1, lambda); 
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/// output parameters into double vector pars used by lowOrg 
kw.rpvDisplay(pars) ; 

/// generate spectrum in RpvSoftsusy object kw 
kw.lowOrgCrpvSugraBcs, mgutGuess, pars, sgnMu, 

tanb, oneset , gaugeUnif ication, ewsbBCscale) ; 

/// outputs for this scan 

cout « lambda « " " « kw.displayPhysO .msnu.displayO) « " # " 
« kw.displayProblemO « endl; 

> 

catch(const string & a) ■[ 
cout « a; exit(-l); 

catch(const char *a) { 
printf(a); exit(-l); 

> 

} 

After including a header file, global variables are defined. These are all described in the Rp manual After 
setting the output accuracy, the program output begins with a title print-out. Then follow some variables speci- 
fying the Standard Model input parameters, the MIXING switch, which determines how any quark mixing is im- 
plemented (as in ref. |5]) and the iteration precision of the output, TOLERANCE. The running masses of the SM 
fermions and the QED and QCD gauge couplings are determined at Mz from data with the method toMz. If the 
switch gaugeUnif ication=true, SDFTSUSY will determine mGutGuess as the scale Mcurof electroweak gauge 
unification. In order to do this, it requires an initial guess, which must be supplied as the initial value of the variable 
mGutGuess (in GeV), and is later over-written by the prograrrQ. 

The next step is the definition of the mSUGRA parameters Ao/GeV=aO, Mi/2/GeV=ml2, tan/3 =tanb and 
mo/GeV=mO. Next, a for loop performs the scan over A'^^-^{Mcut)- In the example given (mSUGRA), the first 
three parameters are pars (1 ) = mO ; pars (2) = ml2 ; pars (3) = aO ; . The pars vector is needed to keep track 
of the boundary conditions set at Mcut- In the iterative SOFTSUSY algorithm the parameters in the RpvSoftsusy 
object change due to the RGE running, they are re-set in every iteration at Mqut from the unchanged DoubleVector 
pars parameters. 

We do not fill the other 102 Ijip entries of pars explicitly. This would be tedious and an additional source of 
potential bugs. Instead, we fill the RpvSoftsusy object itself using the setLsunPrime method in this example. We use 
the rpvDisplay method: this fills the pars vector automatically with what was set already inside the RpvSoftsusy 
object, while leaving the first nine entries in the vector unchanged. The rpvDisplay method automatically changes 
the length of pars appropriately. After this, the actual SOFTSUSY main driving method lowOrg is called, the first 
argument specifying the type of boundary condition (currently rpvSugraBcs), which assumes that pars has already 
been prepared by using the rpvDisplay object. This is followed by the output of the GUT-scale coupling A, the pole 
tau sneutrino mass and any problems in the calculation of the parameter point. Finally, the catch commands print 
any errors produced by the code. 

Appendix D. Object Structure 

We now go on to sketch the objects and their relationship to each other This is necessary information for any 
generalisation beyond the Ijtp MSSM. Only methods and data which are deemed of possible importance for prospective 
users are mentioned here, but there are many others within the code itself. 



If the user wishes to provide this, 2 x 10 GeV is a good initial guess for Mquj. 
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RGE 



QedQcd 




MssmSusy 








RpvSusyPar s 




Sof tParsMssm 








RpvSof tsusy 




MssmSof tsusy 



RpvSoftPars 



Figure D.2: Heuristic high-level object structure of SOFTSUSY. Inheritance is displayed by the direction of the arrows. 



Appendix D.l. Tensor 

The SOFTSUSY program comes with its own linear algebra classes of vectors and matrices (e.g. DoubleMatrix) 
that have been introduced in |5]. New in this version is the class Tensor, given in files tensor . h and tensor . cpp. 
This class has been added to implement the three-index tensors containing some of the Ijip couplings into the program. 
For this reason, the class is specifically designed as a vector of three single objects, each of type DoubleMatrix. 
The dimension of Tensor is (3,3,3). The class also contains linear algebra functions for multiplication, addition 
or subtraction with matrices and vectors via over loaded operators. For more detail we refer the interested reader 
to the technical documentation on the SOFTSUSY web-site. The SOFTSUSY internal representation of the trilinear^^ 
couplings is different to the user interface. While the user interface uses the common Aij^LiLjE e notation, this 
coupling is internally represented in the 1 15] Tensor notation in terms of three matrices as in Eg. IB. 2! Yet, this does 
not stop the user from only operating the code using the user interface and the usual Aijk conventions, /l^.^., /ij^.^ are 
also stored within analogous Tensor representations of three 3x3 matrices. 



Appendix D.2. General object structure 

From an RGE point of view, data in a particular quantum field theory consist of a set of parameters defined 
at some renormalisation scale Q. A set of p functions describes the evolution of the parameters and masses to a 
different scale Q' . This concept is embodied in an abstract RGE object, which contains the methods required to run 
objects of derived classes to different renormalisation scales (their beta functions). The other objects displayed in 
figure ID^ are particular instances of RGE, and therefore inherit from it. QedQcd, MssmSusy, Sof tParsMssm and 
MssmSof tsusy objects encode the Rp part of the MSSM and its SM input data [Ist]. RpvSusyPars contains all of the 
supersymmetric couplings contained within Eq. (|2]l. RpvSoftPars contains the Ijtp soft supersymmetry breaking 
parameters listed in Eqs. Q and (|4]l. RpvSof tsusy is the Ijtp generalisation of the MssmSof tsusy class, and contains 
all Ijip MSSM couplings along with their beta functions. Methods for the RpvSof tsusy class exist to perform the 
minimisation of the neutral scalar potential as well as the calculation of Yukawa couplings described in section |4] 
Neutrino, neutralino and chargino masses and mixings are calculated within this class. Code in the MssmSof tsusy 
class organises and performs the main part of the calculation, using polymorphism to detect the correct p functions 
to use (in this case, the ^^-MSSM /? functions). All of the Rpv objects contain default constructors and destructors, 
as well as overloaded >>,<< operators for input and output. There is always an implicit dependence of running RGE 
quantities on the current renormalisation scale Q. Thus, if a method is called that returns one of the object's couplings 
or masses, that object will return it at the current scale Q of the object. In the following, we provide basic information 
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data variable methods 



Tensor lu, Id, le 


trilinear Ijip 


superpotential 


displayLam, displayLamPrime 


Afi, Ad, Af/ 


couplings 




displayLamPrimePrime 








setLam, setLamPrime 








setLamPrimePrime 


DoubleVector kappa 


bilinear Ijip 


parameter 


displayKappa 


^,=i,2,3/GeV 






setKappa 



Table D. 1: RpvSusyPars class data and accessor methods. See Eq. jB.2t for a translation between the structure of Tensor lu. Id, le. 



data variable 




methods 


DoubleVector mHllsq 
m] /GeV2 


bi-linear scalar 
parameters 


displayMhllSquared 
setMhllSquared 


Tensor her, hdr, hur 

{^,,.,^;,„/^;;,)/Gev 

i,;,^G {1,2,3) 


trilinear .^p scalar 
interactions 


displayHr 
setHr 


DoubleVector dr 

A=i,2,3/GeV2 


bilinear .^p scalar 
parameters 


displayDr 
setDr 



Table D.2: RpvSof tPars class data and accessor methods. 



on the classes associated with Ijip, so that users may program using the class structure of SOFTSUSY. More detailed 
and technical documentation on the program should be obtained from the SOFTSUSY website. 

Appendix D.3. RpvSusyPars class 

Each of the higher level objects described in this appendix have explicitly named display and set methods that 
are used to access or change the data contained within each object. In table ID.lU as in the following tables in this 
section), these accessing methods are listed on the same row as the relevant data variable. The data and input/output 
methods in the RpvSusyPars class are presented in table ID. ll When using Tensor objects lu , Id , le, there exists an 
enumerated type RpvCouplings e {LU, LD, LE) used as arguments to the displayLajnbda,setLajiibda methods. 
This argument selects the type of coupling ^'ij/. or Aijic). 

Appendix D.4. RpvSof tPars cZass 

The data and input/output methods in the RpvSof tPars class are presented in table ID. 21 The displayHr, 
setHr methods take a parameter of the enumerated type RpvCouplings as their first argument to select a partic- 
ular tri-linear scalar interaction that is (h'.'-i^, h'.ji^ or hij^ depending upon the argument). 

Appendix D.5. RpvSof tsusy c/ois 

The data and important methods in the RpvSof tsusy class are presented in table ID. 31 As well as standard 
constructors and destructors, there exists a method beta, that calculates the numerical values of the y6-function. The 
rpvDisplay method is used in the example program rpvmain . cpp, and fills a vector with MSSM running parameters 
in a certain order. rpvSet is used to set MSSM running parameters according to a DoubleVector argument, assuming 
the same order as rpvDisplay. If the user wishes to provide their own function encoding high-scale boundary 
conditions on the soft supersymmetry breaking terms, they must provide a function 

void (*boundaryCondition) (MssmSof tsusy &, const DoubleVector &) 
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data variable methods 



DoubleVector 

v/=i,2,3/GeV 


snuVevs 


sneutrino VEVs 


displaySneutrinoVevs 
setSneutrinoVevs 


DoubleVector 

'Mv,=i,2,3/Ge'V 


nuMasses 


neutrino masses 


setNeutrinoMasses 
displayNeutrinoMasses 


DoubleMatrix 
(7 by 7) 


neutralFermMixing 


neutral fermion 
mixing 


displayNeutralMixing 
setNeutralMixing 


method 




function 




void rpvDisplay (DoubleVector) 


returns all 102 MSSM running parameters as a 
DoubleVector, leaving the first 9 entries untouched 


void rpvSet 


(DoubleVector) 


sets all 102 MSSM parameters from a user-provided 
DoubleVector, Starting from its 10''' entry 



Table D.3: RpvSof tsusy class. 

which is then passed as the first argument to the main driving method lowOrg jst]. One can specify mSUGRA 
conditions for the Rp parts, plus all specified interactions at the high scale by using rpvSugraBcs. Alternatively, 
one can specify 7?^ -conserving gauge mediated supersymmetry breaking conditions plus Ijtp interactions at the scale 
Mmess by using rpvGmsbBcs. In this case, the first elements of pars should contain n, the number of vector like 5-plets 
of messenger fields, the messenger mass scale M„,i,ss in GeV, A and Cgrav, the constant that determines the gravitino 
mass i29ll . rpvAmsbBcs implements the minimal anomaly mediated supersymmetry breaking f3d\ assumption, but 
neglects Ijip couplings in the high-scale boundary condition. For this case, the first two parameters of pars should be 
m32 and mO, respectively. This approximation ought to be reasonable for small dimensionless^^fp couplings. 
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